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Abstract. A systematic study of fractional revival at two sites in XX quantum spin chains 
is presented and analytic models with this phenomenon are exhibited. The generic models have 
two essential parameters and a revival time that does not depend on the length of the chain. 
They are obtained by combining two basic ways of realizing fractional revival in a spin chain each 
bringing one parameter. The first proceeds through isospectral deformations of spin chains with 
perfect state transfer. The second arises from the recurrence coefficients of the para-Krawtchouk 
polynomials with a bi-lattice orthogonality grid. It corresponds to an analytic model previously 
identified that can possess perfect state transfer in addition to fractional revival. 


1. Introduction 

Spin chains with engineered couplings have proved attractive for the purpose of designing 
devices to achieve quantum information tasks such as quantum state transfer or 
entanglement generation m m IE One reason for the interest is that the internal 
dynamics of the chain takes care of the processes with a minimum of external 
intervention required. In this perspective, a desired feature of such chains is that they 
exhibit quantum revival. For perfect state transfer (PST), one wishes to have, for 
example, a one-excitation state localized at the beginning of the chain evolve with unit 
probability, after some time T, into the state with the excitation localized at the end of 
the chain. Such a relocalization of the wave packet is what is referred to as revival M. 
Fractional revival occurs when a number of smaller packets seen as little clones of the 
original one form at certain sites and show local periodicities HIE]. Its realization in a 
spin chain would also allow to transport information with high efficiency via one of the 
clones. Moreover, in a balanced case where there is equal probability of finding clones at 
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the beginning and at the end of a chain, fractional revival would provide a mechanism 
to generate entangled states. It is hence of relevance to determine if FR is feasible in 
spin chains, and if so, in which models. Some studies have established the fact that 
this effect can indeed be observed in spin chains The present paper offers a 

systematic analysis of the circumstances under which fractional revival at two sites in 
quantum spin chains of the XX type will be possible. The main question is to determine 
the Hamiltonians H that will have the fractional revival property. Like for PST, the 
conditions for FR are expressed through requirements on the one-excitation spectrum 
of H. One then deals with an inverse spectral problem that can be solved with the 
assistance of orthogonal polynomial theory. In the full revival or PST case where this 
analysis has been carried out in detail (see [10] for instance), one sees that a necessary 
condition is that the couplings and Zeeman terms must form a three-diagonal matrix 
that is mirror-symmetric. Furthermore, models based on special orthogonal polynomials 
have been found where PST can be exhibited in an exact fashion. These analytic models 
are quite useful; in fact, the simplest one PH was employed to perform an experimental 
quantum simulation with Nuclear Magnetic Resonance techniques of mirror inversion 
in a spin chain ra- We shall here also provide analytic models with FR. 

The outline of the paper is as follows. In Section 2, we present the Hamiltonians 
H for the class of XX spin chains that will be considered. The reader is reminded 
that their one-excitation restrictions J are given by Jacobi matrices and the orthogonal 
polynomials associated to the diagonalization of those J’s are described. In Section 3, 
we review the elements of the characterization of XX spin chains with PST that will 
be essential in our fractional revival study. In particular, the necessary and sufficient 
conditions for PST in terms of the spectrum of J, the mirror symmetry and the 
properties of the associated orthogonal polynomials will be recalled. 

In Section 4, we undertake a systematic analysis of the conditions under which 
fractional revival can occur at two sites. Up to a global phase, the revived states will be 
parametrized in terms of two reals amplitudes sin 9 and cos 9 and of a relative phase ip\ 
the case 6 - 0 will correspond to the PST situation. It will be shown that in general, for 
FR to occur, the one-excitation spectrum of the Hamiltonian must take the form of a 
bi-lattice, i.e. the spectral points need to be the union of two uniform lattices translated 
one with respect to the other by a parameter J depending on 6 and ip. Two special cases, 
namely if - 0 and ip = n/2, will be the object of the subsequent two sections. In the 
first case the spectrum condition is the PST one and in the second case, it is rather the 
mirror symmetry that is preserved. These two cases will provide the ingredients of a 
two-step process for obtaining the most general XX chains with FR at two sites. 

In Section 5, it will be shown how spin chains with FR can be obtained from 
isospectral deformations of spin chains with PST. This has the relative phase ip = 0. 
Analytic models with FR will thus be obtained from analytic models with PST. It will 
be observed that only the central parameters of the chain will need to be modified 
so as to make FR happen. Mirror-symmetry will be seen to be replaced by a more 
complicated inversion. The orthogonal polynomials corresponding to the deformed 


Quantum spin chains with fractional revival 


3 


Jacobi matrix will be shown to have a simple expression in terms of the unperturbed 
orthogonal polynomials associated to the parent PST Hamiltonian. Their knowledge 
will be relevant to the construction, at the end of Section 6, of the FR Hamiltonian with 
arbitrary parameters 9 and 

Section 6 will deal with the special case when the relative phase is 7r/2, that is when 
one amplitude is real and the other purely imaginary. It will be seen that J is mirror 
symmetric in this situation with its spectrum a bi-lattice. As shall be explained, an exact 
solution to the corresponding inverse spectral problem turns out to be already available 
and is provided by the recurrence coefficients of the para-Krawtchouk polynomials 
introduced in [13j. It is remarkable that these somewhat exotic functions naturally arise 
in the FR analysis. We shall demonstrate that for certain values of the parameters, 
the models thus engineered possess both FR and PST. We shall finally come to the 
determination of the generic Hamiltonian with FR and show that is is obtained from the 
recurrence coefficients of para-Krawtchouk polynomials perturbed in the way described 
in Section 5; in other words it is constructed by performing an isospectral deformation of 
the Jacobi matrix of the para-Krawtchouk polynomials. The two arbitrary parameters 
are related to the the bi-lattice and to the deformation parameter. We shall summarize 
the outcome of the analysis and offer final remarks to conclude. 


2. XX quantum spin chain models with non-uniform nearest-neighbor 
interactions 


We shall consider XX spin chains with N + 1 sites and nearest-neighbor interactions 
that are governed by Hamiltonians H of the form 


N-l 


N 


H = 




i U U1 


v v 

+ 


) ++ i) 


e=o 


£=0 


( 2 . 1 ) 


that act on (C 2 )® iV+1 . The Ji are the coupling constants between the sites £-1 and £ 
and the Bi are the strengths of the magnetic fields at the sites £, where £ = 0,1,..., TV. 
The index £ on the Pauli matrices a x , a y , a z indicates on which of the (AW 1) C 2 factor 
these matrices act. If | t) and | |) denote the eigenstates of a z with eigenvalues +1 
and -1 respectively, a x and a y are known to act as follows in that basis: a x \ f) =| I), 
a x | |) =| |); cr y | t) = ~i | I), <J V | I) = i\ t)- The Hamiltonians H are invariant under 
rotations about the z-axis, i.e. 


1 N 

[if,-J> ( = + l)]=0, (2.2) 

z e=o 

and as a consequence, the eigenstates of H split in subspaces labeled by the number of 
spins over the chain that are in the state | t), which is a conserved quantity. In the 
following we shall focus on the subspace with one excitation and we will denote by J 
the restriction of H to that subspace, equivalent to C N+1 . A natural orthonormal basis 
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1^) - (o, o,..., 1,..., 0) T , e = 0,l,...,N, (2.3) 

where the only “1” in the entry corresponds to the only state | t) being at the site t. 
The action of J on these basis vectors is directly obtained from (12.11) and given by 

J\£) = J M \£ + 1) + B e \£) + Jt\£-1), (2.4) 


where it is assumed that Jo = Jn+i = 0. The restricted Hamiltonian J thus takes the 
form of a (N + 1) x (N + 1) three-diagonal Jacobi matrix 


J 


( B 0 Ji 
J\ Bi J 2 

J2 B2 

Jn 

{ J N B N) 


(2.5) 


It is known that Jacobi matrices are diagonalized by orthogonal polynomials. Let us 
record, in this connection, results that will prove useful in our study. Let |A) be the 
eigenvectors of J: 


r/|A) — A | A). 


( 2 . 6 ) 


Since J is Hermitian the eigenvalues A are real. We shall assume that they are non¬ 
degenerate, i.e. that they take N + 1 different values A s for s - 0,1,..., N. This is the 
case when the entries Jp are positive. We will take the eigenvalues to be ordered 


A 0 < Ai < ••• < A 


N • 


Consider the expansion of |A S ) in terms of the basis vectors \d), 

N 

lA.HEWrflQ 

e=o 

and write the elements W s g of the transition matrix in the form 


(2.7) 


( 2 . 8 ) 


w B t = ir.n x.-( .V.) = \ZwIxf(A s ). 


It follows from 112 fit that Xr(A s ) are polynomials satisfying the three-term recurrence 
relation 

Je+i Xe+i(^s) + Bp X£(A s ) + JeXe-i(^s) = Ax^(A s ), (2.9) 

with initial condition 


Xo = 1, 


X-i = o. 


( 2 . 10 ) 
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Since both the eigenbasis {|A S )}^ 0 and the basis {|£)}^ 0 are orthonormal and all the 
coefficients are real, the matrix (W) s t is orthogonal and the condition {£\P) = 
implies that the polynomials \y(A) are orthogonal on the finite set of spectral points Ay 


N 

^ a sXm(A s )y ri (A s ) — S mn . 

s=0 

Note that the weights w s are normalized, that is 

N 

E w s = !, 

s =0 


( 2 . 11 ) 


( 2 . 12 ) 


since Xo(A s ) = 1. Owing to the orthogonality of W, in addition to (12.8ft . we also have 


N 


N 


|£) = Z W si \K) = Y,V^sXn(\s)\K)- 


s =0 


s=0 


It is sometimes useful to use the monic polynomials 
Pe(X) = yheXe(X), yhe = JiJr-'Je-, 


(2.13) 


(2.14) 


whose leading coefficient is 1: Pe(\) - X e + 
characteristic polynomial 

Pat +1 (A) = (A-Ao)(A-A 1 )-(A-A A t). 


The spectrum of J is encoded in the 


(2.15) 


The following formula from the standard theory of orthogonal polynomials gives the 
weights w s in terms of P^( A) and Pn+i(X)' 

hN 


w, = 


s = 0,1 




PN(K)P' Ntl ( A,)’ 

with P( r+1 (A) denoting the derivative of P/v + i(A) with respect to A [14]. 


(2.16) 


3. A review of perfect state transfer in a XX spin chain 

Perfect state transfer is achieved along a chain if there is a time T for which 

e'* TJ |0) = e i4> \N), (3.1) 

where is a real number. In this case the initial state with one spin up at the zeroth 
site will be found with unit probability after time T in the state with one spin up at the 
site N. One then says that the state | t) has been perfectly transferred from one end of 
the chain to the other. The requirements for PST have been determined from studying 
the implications of (13.1ft . Since this will be the backdrop for our discussion of fractional 
revival, it is pertinent to summarize the PST analysis here. Using the expansion (12.13[1 
on both sides of (13.1ft gives 


e- l<p e- lTX > = xn(Xs) 


s - 0 , 1 ,..., N. 


( 3 . 2 ) 
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Xn( X s ) = ±1. (3.3) 

Now because the zeros of Xn(^) must lie between those of xat+i(A) which are located at 
the A s , we must conclude that Xn(A s ) will alternate between +1 and -1 [13]. Given that 
the eigenvalues are ordered, the condition that the weights, as given by (12.161) . must be 
positive leads finally to 

X n(A s ) = (-1) n+s , s = 0,1,... ,N. (3.4) 

In view of (13.41) . we see that (13.21) imposes the following requirement on the spectrum 
of J: 


0 -jTA» 


= e^(-l) 


N+s 


s = 0,1,... ,7V, 


(3.5) 


which is tantamount to the condition that the successive eigenvalues are such that 

A s+ i - A s = —M s , (3.6) 

with M s an arbitrary positive odd integer. We have thus with (13.41) and (13.5ft . the 
necessary and sufficient conditions for PST to occur. It is seen from (12.16j) that the 
necessary condition (13.41) is equivalent to requiring that the polynomials associated to 
the Jacobi matrix J be orthogonal with respect to the weights 




Pn+ i(As) 


s = 0,1,..., TV. 


Remarkably, the weights given by (13.71) have the general property [15] 


Y, W2 * - Yj W2s + i - 9 - 

S S " 


(3.7) 


(3.8) 


In turn, (13.41) or (13.61) can be shown to hold if and only if the matrix J is mirror 
symmetric with respect to the anti-diagonal [IQ], that is if and only if 


RJR = J, 

with 

n 

R= 1 

a t 


(3.9) 


(3.10) 


In terms of the couplings and magnetic field strengths, one sees from ( 12.51) that this 
symmetry amounts to the relations 


Jn ~ Jn+I-ui 


Bn ~ 


(3.11) 
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Matrices with that property are referred to as persymmetric matrices. A simple and 
direct proof that (13.11) implies (13.91) can be found in [2]. 

Let us now observe that this required symmetry will lead not only to PST but to 
a complete mirror inversion of the register at time T if the spectral condition (13.5ft is 
satisfied. First note that R is an involution, that is R 2 - 1. Because J is reflection- 
invariant, its eigenstates can be taken to also be eigenstates of R and since we have 
assumed that the spectrum of J is not degenerate, each eigenstate must hence be of 
definite parity e s , equal to either +1 or -1. We hence have 

R\X s ) = e s \X s ). (3.12) 

From (12.81) . we have 

N N 

R 4) = E \Ak Xi(X s ) \N-£) = Y J y/wiX n-i(Xs) 14 (3.13) 

£=o e=o 

but given (13.121) . using (12.81) again we find that 

XN-e(X s ) = ( -l) N+s xt(X s ), (3-14) 

where we have determined e s by setting l - 0 and using (13.41) . As a result, we see that 

N N 

(k\ e~ lTJ \£) = E e~ lTXs w s xe(X s )xk(X s ) = e l<t> Y i (-l) N+s w s xe(X s )xk(X s ) 
e=o s=o 

N 

= e l4> E XN-i(Xs)xk(X s ) = e^S^N-e, (3.15) 

s=0 

with the successive help of (j3.5[) . (13.141) and (12.111) . In other words, there is probability 
1 of finding at the site N - n, after time T, the spin up initially at site n. In matrix 
form, we have shown that (13.4[) or (13.6[) and (13. 5 p have for consequence that 

e~ iTJ = e^R. (3.16) 

Now the main issue is to determine the Hamiltonians H of type (12.11) that have the 
PST property. This is an inverse spectral problem since we start from conditions on 
the eigenvalues. As it turns out, this specific problem has been well studied [16]. 
The outcome is that once a spectrum satisfying (13.5[) or (13.61) is given as input, the 
Hamiltonian with the desired properties is uniquely determined. This stems from the 
fact that the weights w s entailing mirror-symmetry are uniquely prescribed by (13.71) 
and that the corresponding orthogonal polynomials can be unambiguously constructed. 
Their recurrence coefficients henceforth provide the couplings J\ and the local magnetic 
fields If of a Hamiltonian H with PST. One algorithm for obtaining the orthogonal 
polynomials is described in mg. By considering natural types of spectrum, analytic 
models of spin chains with PST have thus been obtained and associated to various 
families of orthogonal polynomials DU: Krawtchouk and Dual Hahn HU, Dual -f 
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Hahn using lag, g-Krawtchouk |2f and (/-Racah ra- Of particular interest for what 
follows are models connected to the so-called para-Krawtchouk polynomials pT3j . The 
paradigm example of analytic models, investigated in m, is obtained by considering 
the linear spectrum 
7 r 

t 1 ' 

for which (13.71) yields the binomial distribution 
TV! / i' JV 


A, = 7_(s-JV/2), s = 0,l,...,JV, 


(3.17) 


= 


ii(s) 


(3.18) 


s\(N -s)\ 

The associated polynomials are known to be the symmetric Krawtchouk polynomials 
with recurrence coefficients 


Bp - 0, 


2 _ 7r 2 £(N + l-£) 

'h ~ 


p 4 • ( 3 - 19 > 

The mirror symmetry (13.111) is manifest. This is one instance where the one-excitation 
spectrum dynamics is exactly solvable. The general transition amplitude from state \£) 
to state | k) in time t has been calculated in [22] to be 


{k\e 


-itj 




-i — t\k+£ 


r 


(l + e" 1 -*) 


-i—t\N-k-£ 



-4e 


-i — t 


(3.20) 


( l - e -^)2 

where 2 Fi is the classical hypergeometric series [23] . 

It must be stressed hnally that a myriad of analytic models can be generated from 
those associated to known orthogonal polynomials by a procedure known as “spectral 
surgery” [10]. Indeed, it has been shown that the first or last eigenvalues (Ao or Aa r) 
or a pair of neighboring spectral points (A*, Aj+i) can be removed without affecting the 
PST properties [TO]- If Bp (A) are the monic polynomials associated to the original chain 
with Hamiltonian J, the removal of one level, say A*, will give a new Jacobi matrix J’ 
with entries given by 




72 

J £ 1 


Bp - Bp + 1 + A 


h+ 1 


■ Ap 


where 


Ap - 


Bp+\ (A,) 


P t ( A*) 

The monic polynomials corresponding to J' are 
P £+1 (x) - ApPp(x) 


P'M - 


X - Xi 


(3.21) 


(3.22) 


(3.23) 


Formulas (I3.2ip and (j3.23jl can be applied iteratively to obtain new analytic models 
from known ones. The removal of pairs of neighboring spectral points is required in the 
bulk to ensure the positivity of the weights. 
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We are now ready to examine the requirements for fractional revival in a spin chain of 
type XX. We shall limit ourselves to revivals occurring at only two sites chosen to be 
£ - 0 and £ - N, the two ends of the chain. We shall here analyze systematically the 
conditions that will ensure that for a time T 


e~ iTJ |0) = ct |0) + /31 AT), (4.1) 

with \a\ 2 + |/3| 2 = 1. When this is so, the spin up at site £ = 0 at t = 0 is revived at the 
sites £ = 0 and £ - N after time T with amplitude a and /3, respectively. Let us write 

a - eX' sin 26, ft = e^W) cos 26, (4.2) 

with e l< £ and e*h the global and relative phase factors of the two complex numbers a and 
f3. Taking < 6 < | and 0 < < n with </> e R covers all possible amplitudes. The use 

of (12.13ft leads to the relation 

e~ iTXs _ e ict> [ gin 20 + e # cos 26 Xv(A s )], s = 0,1,..., N. (4.3) 

This implies that sin 26 + e*^cos2 &Xn(^s) has modulus 1, that is 

Xn(^s) + 2 tan26* cos-0 - 1 = 0. (4.4) 


Obviously when 6 - 0, we recover condition (13.31) . Equation il 1. Hi indicates that Xat(A s ) 
will take one of two values, 7 and -4, with 7 satisfying 

7 -= -2 tan 26 cos fj. (4-5) 

7 

By a continuity argument as 6 0, upon comparing with (j3.4[) and assuming that 7 is 

the positive root (there must be one because of the interlacing property of the zero of 
orthogonal polynomials), one concludes that, for N odd 


Xn(^2s) 



Xn(^2s+i) - 7) 


for N even 

Xn( A 2s ) = 7 , 


Xv(A2 S+ i) 


1 

7 


(4.6a) 


(4.6b) 


for all s e {0,..., N}. Assume for now that the relative phase factor eis generic. Once 
(14.51) and (14.61) are satisfied, it must still be ensured that (14.31) is obeyed. When N is 
odd, this amounts to the conditions 


cos(TA 2s + </>) = sin 26 -cos 2 d cos?/;, sin(TA 9 S + </>) = — cos 26 sin?/;, 

7 ^ 7 

cos(TA 2s+ i + </>) = sin 26 + 7 cos 26 cos 'if, sin(TA 2s+ i +</>) = -7 cos 26 sin if. 


(4.7a) 

(4.7b) 
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Define £,77 e [0,27r] by setting the right-hand sides of equations (I4.7a.p . (I4.7bl) to be 
respectively cos£, sin£, cos?/ and sin 77 . This implies 

{TX 2s + cf} £ {£,£±27r,£±47r,...}, {TX 2s+1 + (f} c { 77,77 ± 2 tt, 77 ± 4vr,...}. (4.8) 


When N is even, it is seen that the roles of £ and rj are interchanged. On thus observes 
that the spectrum of J must be the following bi-lattices: 


= E + s Q (S - 1)(1 -(-!)*), 

7T 7T 2 

for s = 0,1,... ,7V with 

f£ for N odd 

P = \ , 

1 77 for N even 

and 

f2 + ^(77 — £) for N odd 
I ^(£ - 77 ) for N even 


(4.9) 


(4.10) 


(4.11) 


We have used the latitude in picking the initial point in the sets (14.81) and those obtained 
by the exchange £ <->• 77 for N even, to ensure that 6 is positive. Let us point out that 
the spectra can in fact also be subsets of the bi-lattices (14. 9 1) after appropriate surgery. 

There are two special cases for the phase if that are of particular interest and that 
will be the object of the next sections. Each preserve one of the two necessary and 
sufficient conditions (13.41) and (13.51) for PST. These two distinguished cases are if - 0 
and if = 7 r/ 2 . Let us make initial observations about what happens for those values. 

1 -^ = 1 

It is readily seen from (14.51) and (14.61) that 7=1 and that Xn(X s ) = (-l) iV+s . Hence 
mirror-symmetry is maintained in this case. This is the case considered in [ 8 ]. Equations 
(14. 71) give that 

£ = - 77 = |-20. (4.12) 

The spectral points thus form bi-lattices of the form (I4.9[) with 

4 9 

5= 1±—. (4.13) 

7 r 

The upper/lower sign in the above equation corresponds to N being odd/even. Note 
that 5 e [0,2] for -7 < 9 < when 9 = 0, it is the PST situation, 5 = 1 and (14. 9 p 
becomes the linear spectrum of the Krawtchouk polynomials with (f - ± 1). As 

shall be explained in Section 6 , the model corresponding to the spectral conditions for 
if = | is analytic and can exhibit both FR. and PST. 
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ii. if - 0 

In this case (14. 5 p gives for 7 

7 = cot (^ ~ $) j (4-14) 

and in view of (14.61) mirror symmetry is broken. However using (14.14j) and (14.61) it is 
immediate to check that (j4.3[) becomes 

e ~ iTX ° = e i *(-l) N+ f (4.15) 

which is the spectral condition (I3.5[) for PST. So when if - 0, PST is absent for lack 
of mirror symmetry but the spectrum remains unchanged. This isospectral situation is 
discussed next. 

5. Isospectral deformations of chains with perfect state transfer 

We shall now describe the spin chains with fractional revival that can be obtained by 
isospectral deformations of chains with PST [9], This picture arises when the relative 
phase if is nil. It should be stressed from the outset that the procedure will generate 
analytic models with FR when it is applied to spin chains for which PST can be exactly 
demonstrated. In this section, for the sake of clarity, we shall denote by J the one- 
excitation Hamiltonians of spin chains with FR and by J those of spin chains with PST. 
When if - 0, the FR condition (14.11) reads 

e -* Tj |0) - e i(t> [sin 29 |0) + cos 29 |7V)] (5.1) 

and (14.31) becomes 

e -i<t> e -iT\ s _ s j n 20 + cos 29 xn(^s)- (5.2) 

We observed that because the right-hand side of (15.21) is real, the condition on the 
spectrum of J is the same as for PST, that is (13.51) . For one such spectrum, still 
assumed to be non-degenerate, it must be possible to relate by a conjugation the Jacobi 
matrix with PST to the one with FR. Recall that the PST matrix J can be uniquely 
constructed from the data. There is thus an orthogonal matrix U such that 

J = UJU J . (5.3) 

It then follows that 

e~ iTj = Ue~ iTJ U T = e^URU 1 = e^Q. (5.4) 

Note that the action of Q on |0) is prescribed by (15.11) . This similarity transformation 
is easily found and can be presented as follows. For convenience write U in the form 


U = VR. 


(5.5) 
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/ sin 6 


cos 6 \ 


sin 6 cos 9 
cos 9 - sin 9 


^cos 9 

and for N even, let 
f sin 9 


- sin 9 j 


cos 9 


sin 9 

0 

cos 9 

0 

1 

0 

cos 9 

0 

- sin 9 


^cos 9 


- sin 9) 


(5.6) 


(5.7) 


It is immediate to check that V - V T and that V 2 = 1. It then follows that UU T = 1. 
Note also that detR = -1. Obviously V(0) = R. The matrix Q introduced in (15.4p is 
thus given by 


VRV = Q, (5.8) 

and is obtained from V by substituting 9 by 29 in (15.61) for N odd and in (15.71) for N 
even. Obviously Q 2 = 1. Recall that the PST matrix J is persymmetric RJR = J. It is 
then easy to see that for 

J=UJU T = YJ Vl (5.9) 


condition (15.11) is satisfied. In fact, not only is this realized but in view of (15.41) and the 
expression for Q , we shall have fractional revival between the mirror-symmetric sites t 
and N - £ since we have 


N odd N even 


e- iTJ = 


" sin 29 \£) + cos 29 \N - 1) 
- sin 29 \£) + cos 29 |IV - l) 


£<^ 




and for N even 


-iTJ _ e i0 |D 


(5.10) 


(5.11) 


To sum up, we have seen that a chain with FR can be obtained from any chain with 
PST by conjugating the Jacobi matrix of the latter according to (15.9[) . The resulting 
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operator J is not mirror-symmetric but is seen to be invariant under the one-parameter 
involution Q, that is 

QJQ = 7 . (5.12) 


It is remarkable that the only modifications or perturbations in the couplings and 
magnetic fields that arise when passing from J to J occur in the middle of the chain. 
Indeed, upon performing the conjugation (15.911 with (15.61) or (15. 7 j) . recalling that J is 
persymmetric, one finds that the only entries of J that differ from those of J are 


Jn+i = J n +i cos 29, 

~ 2 2 . (5.13a) 

B jvti = B n -i ± J n +i sin 2d, 

2 2 2 

for N odd and 

Jn = Jn_ (cosd + sind), 

_ 2 2 (5.13b) 

Jn ., = Jn (cos d - sin d), 

2 2 

for N even. When N is even, only the couplings between the three middle neighbors 
are altered. When N is odd, it is only the coupling between the two middle neighbors 
that is affected together with the magnetic held strengths at those two middle sites. 
Note that if all the B£ of J are initially zero, J will only have two Zeeman terms of 
equal magnitude and opposite sign at £ = ^rr- and £ = The fact that in these 

models so few couplings or held strengths of the PST chain need to be adjusted to 
obtain the chain with FR could prove to be a practical advantage. One can imagine 
that the calibration would hrst be done by engineering the couplings so as to reproduce 
the PST mirror inversion and that thereafter the transformation to the FR mode would 
not be technically too prohibitive. It is also interesting to remark that it is possible to 
have no (zero) coupling between two equal parts of the chain and hence two separate 
chains in fact, and yet to keep some transport. Indeed, when 9 - for N even we have 
J n +i = 0. Bn±i = Bn^i ± J n +i and when N is odd, Jn = 0. Jn..-, = V2Jn_. It is clear that 
analytic spin chain models with fractional revival can be obtained from the analytic 
models with PST that are known by performing the isospectral deformations that we 
have described in this section. Take again for example the system associated to the 
Krawtchouk polynomials. Starting with the couplings J^ and magnetic held Bp given in 
(13.19[) and modifying them according to (15.13|) will yield a rather simple Hamiltonian H 
(with one-excitation sector J) with fractional revival. The exact solvability properties 
of the perturbed model will be inherited from those of the Krawtchouk chains. For 
instance, the general transition amplitude between the one-excitation states |£) and | k) 
during time t under the evolution governed by J, that is (A:| e~ lt ^ \l) can be obtained 
directly from the corresponding quantity associated to J and given in (13.201) . Indeed, 


(k\ e- itJ \£) = (k\ Ve~ iTJ V \£) = £ V mk V n , (m\ e~ itJ \n), (5.14) 

mn 
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which will yield a sum of (at most) four terms owing to the special form of V. Let 
us mention that some of the couplings (15.13ft have appeared in studies of entanglement 
generation. The cas 6 = n/8 in [21] and the case N even in [2]. 

To complete the discussion, we shall conclude this section by providing information 
on the relation that the orthogonal polynomials associated to the Jacobi matrix J = V JV 
with fractional revival bear with those attached to the matrix J with PST. This will 
offer consistency checks and will be of relevance when considering the generic situation 
when the relative phase ^ of (14,31) is arbitrary. Let |A S ) be the eigenstates of J 

J\\ s ) = \ s \\ s ). (5.15) 


Recall that J and J have the same spectrum. We have an expansion analogous to (12. 8 1) 
in terms of a different set of orthogonal polynomials yy(A): 


~ N 

|A a ) = Y,>/filxiM\Z), 

£=0 

where the weights w s are given by the formula (12.161) that now reads 

~ h N 

w s = — -—-, 

Pn(\s)P^ +1 M 


(5.16) 


(5.17) 


with the monic polynomials P\ dehned by 


P, 




JiJt-Jz- 


From (15.131) we see that 


(5.18) 


\ 


hN 

-— = cos 26/. 
llN 


(5.19) 


Now, in view of (14. 6 j) and (13. 7 j) , we find that the weights iv s and w s are related as follows: 
| 7 cos 29 W‘ 2 S N odd 


W 2s 


- cos 29 1 U 2 s N even 


W2s+1 


A cos 29 w 2 s+i N odd 
I 7 cos 20 w 2 s+i N even 


(5.20) 


Since 7 + 7 1 = 2 sec 29 as is readily observed from (14.140 . one checks in particular that 
= 1 using (I3.8j) . With |A S ) = R|A S ), using (12.81) . one also has 


N 


IA s ) = V |A S ) = Y, V™lXk(K)Vik 1^), 

£,k =0 

and upon comparing (15.161) with (15 .21 p . one hnds the relation 

N 

Vti~sXe(K) = Y Y^VikXki As). 


(5.21) 


(5.22) 





















Quantum spin chains with fractional revival 15 

Given the form of V and the property (13.141) of the polynomials, one can write 

VwsXi(^s) = y/uhiVtj + (-l) N+s V^N-e)xe(K)i (5.23) 

together with 

V^sXK (K) = Vw~sX*( K), (5.24) 

for N even. In this last instance, note that (13.141) gives 

X*(A fl ) = (-l) s x*(A a ), (5.25) 

meaning that Xiv(A) is zero on all the odd eigenvalues: Xjv(A 2s+ i) = 0. The simple 
2 2 
trigonometric identities 


7±7 = 


/ 7T . 

A / 7 T 

\ (2 sec 29 

[--6 

1 ± cot ( — 9 

= ^ 

V 4 

) \ 4 

/ (^-2 tan 26 


(5.26) 


allow to show that 

7 cos 29 - (sin 9 - cos 9) 2 y _1 cos26> = (sin# + cos(9) 2 (5-27) 

Using these relations and examining each case, one checks that (I5.23j) and (15.24[) imply 

«(A.) = xe(\.), e < 

iV + 1 (5.28a) 


Xr(A 2s ) = 7 1 X4(A 2s ), 
X^(A 2 s+ i) = 7X^(A 2s+ i), 
for N odd and 

Xi{ As) = Xi(^s), 

_ X4(A 2s ) 

X4.(A 2s ) = 


9 

1V+ 1 


N 

< 2 ’ 


sin 9 + cos 9 ’ 
Xjv(A 2s+ i) = Xjv(A 2s+ i) = 0, 


t -P 

2 


N 

~2' 


(5.28b) 


_ N 

xi(^2s) = 7X£(a 2s ), 

X<(A 2a+ i) = 7 _1 xKA 2s+ i), t > y,, 

for N even. The proper sign should be chosen in taking square roots of the relations 
(15.271) so that (14.61) is fulfilled. 

Consider now the difference between the monic polynomials Pjj and P^: 


Pn - Pn - \/^iv(cos 29 \ N - Xn)- 


(5.29) 
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Pn(^ 2 s) - ^at(A 2s ) = V^v(l -7 1 cos20), (5.30) 

Pn(^ 2 s+i) - Pn(^ 2 s+i) = \fhN (7 cos ‘Id - 1 ), (5.31) 

for N odd. The right-hand sides are interchanged for N even. It is readily checked with 
(15.261) that 

(1 - 7 -1 cos 20) = 7 cos 20 - 1. (5.32) 

Hence for N odd and even, 

Pn(K) - PnM = \//pv(7cos 26 - 1 ). (5.33) 

This shows that Pjv ~ Pn, a polynomial of degree N - 1, is equal to a constant for N 
points and must hence be identically equal to that constant. We thus have 

-f)v(A s ) = Pjv(A s ) + Coj (5.34) 

with 


Co = cos29 - 1 ). (5.35) 

From the knowledge of Pat + i(A s ) = P/v+i(A s ) and of Pjv(A s ), using the recurrence 
relations, it is possible to show by induction that 

Pi - Pi, (5.36a) 

and 


Pn-i - Pn-i + Ce.Pi, 


(5.36b) 


for 

f0,...,^± N odd 
(0,..., 4 - 1 N even 

with in addition 


(5.37) 


^V(A) = Pn(\) = (A - A X )(A - A 3 )---(A - Ajv-i), 
when N is even. The constant Q in (I5.36b|) is given by 
Co 


U 


J 2 ... /2 /2 ' 

J N+1-1 J N-1 J N 


(5.38) 


(5.39) 


Details will be given elsewhere |T5j. The fact that the polynomials Pi are equal to 
the unperturbed polynomials Pi for the first half of the indices/degrees was expected 
because the recurrence coefficients are the same up that point. It is readily checked that 
the evaluations (I5.28P on the spectral points are entirely consistent with the formulas 
(15.361) when one allows for (I2.14|) . 













Quantum spin chains with fractional revival 

6. The bi-lattices models and para-Krawtchouk polynomials 


17 


We saw in Section 4 that within the class of XX spin chains with non-uniform nearest 
neighbor couplings, the general conditions in order to have fractional revival at two 
sites are two-fold. One, the spectrum {A s } of the one-excitation Hamiltonian J must 
be comprised of the points of the bi-lattice (14.911 or of an ordered subset of those grid 
points resulting from the removal of consecutive eigenvalues. Second, the transition 
matrix that diagonalizes J must be made out of polynomials that are orthogonal with 
respect to the weight w s given by 




U>2s = 


_ iVhN 

p n + i(*2.) 

tPW(^) 


N odd 
N even 




W2s+1 = 


Vhisr 

T' P at+i( A 2s+1 ) 
7 V^N 

p n+ i(-Wh) 


N odd 
N even 


( 6 . 1 ) 


with 7 the positive root of (14. 5 1) and P' N+l (\) as before, the derivative of the 
characteristic polynomial. Finding the specifications of the corresponding spin chain 
amounts to an inverse spectral problem that can be solved by finding the polynomials 
occurring in the transition matrix and thereafter their recurrence coefficients which are 
the entries of J. In Section 4 still, we pointed out that there is an interesting special 
case that arises when the relative phase if - |. When this is so, 7 = 1, w s is given by 
(13.711 and we know that J is mirror-symmetric. This case has been considered in [ 8 ] 
and we shall discuss it in detail here. The authors of [ 8 ] have determined numerically 
the persymmetric Jacobi matrix in the perfectly balanced situation 9 = tc/8. We shall 
indicate that there is in fact an exact description for any 9. Indeed, chains with bi¬ 
lattice spectra and mirror-symmetric couplings have been analyzed with the help of 
the para-Krawtchouk polynomials that two of us have identified and characterized in 
|l~3j . Since their Jacobi matrix is persymmetric, these models are poised to admit 
PST. The circumstances under which they shall exhibit PST in addition to FR will 
be discussed. We shall conclude the section by returning to the general case. We 
shall explain that it can be realized by combining the construction of the persymmetric 
matrices associated to bi-lattices and para-Krawtchouk polynomials with the isospectral 
deformations described in Section 5 and possibly surgeries. 

When if = 7t/ 2, the spectral condition (14.31) becomes 


e -iT\, =e it [ sin 20 + (—!)»«; cos 20 }, 


( 6 . 2 ) 


which amounts to 


e -iT\2s = e i<j> e i( f-20), e -iT\ 2a+ i = e # e -i(f-20)_ 


(6.3) 


Let us now calculate the transition amplitude (k\ e lTJ \l) in analogy with what was 
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(k\ e~ iTJ \£) 

= E e_ * TA2S W 2sXi{^2s)Xk{Ms) + Y e~ lTX2a+1 W2s+lXe(^2s+l)Xk(^2a+l) 


2s 


2s+l 


. e ^ e ™l 2 \cOs26(Y W 2sX(.{^2s)Xk{^2s) ~ Y w 2s+lXi(^2s+l)Xk(>^2s+l)) 

L V 2 s 2s+l ' 

i sin 29 (l W2sXl{^2s)Xk{Ms) + X! W 2s+lXe(^2s+l)Xk(^2s+l)^ 


(6.4) 


2s+l 


Making use of (13.141) . one finds that 

(k\ er lTJ \£) = e^ sin 26 + i cos 268N-e,k] , (6-5) 

which shows that 

e~ iTJ \£) = [sin 26\d.) + i cos 26 \N - £)\ , £ = 0,1,..., Ah (6.6) 

As observed in [8], we see that a state localized at site £ will be revived at the sites £ 
and N - £. In matrix form, we have found that 


e -iT J 


/ sin 26 

i cos 2 fh 

\i cos 26 

sin 26 j 


(6.7) 


for N odd. In the same way, one shows that for N even 
f sin 26 i cos 26^ 


e -iTJ _ e i(j> 


e i(f-20) 


( 6 . 8 ) 


^zcos 26 


sin 26 t 


The set of couplings and magnetic held strengths for which Hamiltonians of the form 
(12.11) will lead to this behavior has been provided explicitly in [TS]. They happen to be 
formed of the recurrence coefficients of the orthogonal polynomials that have been called 
the para-Krawtchouk polynomials. These OPs are precisely those that are associated 
to persymmetric Jacobi matrices with the bi-lattice spectra 

x, = s + i(i-l)(l-(-l)*). (6.9) 

They have been constructed with the help of the Euclidean algorithm, described in [10], 
from the knowledge of the two polynomials Pn+i and Pn, the former being prescribed 
by the spectrum and the latter by the mirror symmetry. They have been named 
para-Krawtchouk polynomials on the one hand because their spectrum coincides, when 
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N -*■ oo, with that of the parabosonic oscillator [25] and on the other hand because they 
become the standard Krawtchouk polynomials when 5=1. Their recurrence coefficients 
and properties are given in ra One has for N odd 


N -1 + 8 

Jjf - -z-, 


1 


A 


£(N+l-£)((N+l-2£) 2 -8 2 ) 
(N -2£)(N -2£ + 2) : 


and for N even 

- N-l + 6 . 0-l)(JV + l) 

■tie = ---^-:- 


1 


, 21 - N - 1 21+1-N 


)■ 


Je = 1 


A 


£{N + 1 - £)((2£ - N - l) 2 - (5 - l) 2 ) 
(2£ - N - l) 2 


( 6 . 10 ) 


( 6 . 11 ) 


for £ - 0,1,..., N. These formulas correspond to the lattice (16.91) . It is easy to see from 
the recurrence relation 


xP e (x) = P e+ i(x) + B e P e (x) + f t Pi- i(x), 


( 6 . 12 ) 


of the monic polynomials for instance, that an affine transformation of the lattice points 


x s - ax s + b , 


(6.13) 


will lead to orthogonal polynomials Pg(x) with recurrence coefficients given by 


Bp — aBg + 5 , Ji — dJ t- 


(6.14) 


Note that the diagonal terms Bp, that is the magnetic fields, are the same at every site 
for N odd; see (16.101) . They can thus be made equal to zero by an affine transformation. 
This is not so for N even however. Comparing (I4.9|) with (16.91) and making use of (I6.14[) . 
we see that by choosing the global phase f to be 


f>=\{N-l + 5) 


N odd 
^~ 1 ' ) N even 


(6.15) 


and in view of (14.!()[) - (14.13[) . the spin chains with the fractional revival features described 
in this section have their couplings and magnetic fields given by 

Be- 0, Je-—Ji , (6.16) 

for N odd and 


B,= 


e_ 

T 


(JV + 1)( 


1 


2£ - N - 1 2£+l- N 


y 


j ,=i j t , 


(6.17) 


for N even, with given by (16.101) and (16.111) where 5 = 1 + AO/tt for N odd and 
5=1- 4:6/n for N even. As observed also numerically in [8], relative to the coefficients 































Quantum spin chains with fractional revival 


20 


of the Krawtchouk chain given in (13.191) . the magnetic fields remains zero for N odd 
while they are proportional to 9 for N even. Note that the Krawtchouk chain parameters 
are recovered when 9-0 and that contrary to the isospectral models with fractional 
revival covered in the last section, here, all the Ji are modified in comparison with those 
of (13.191) . Interestingly, the para-Krawtchouk models have been shown in [13] to enact 
PST for 


M 1 


(6.18) 


where Mi and M 2 are positive co-prime integers and M x is odd. Let us here explain 
how spin chains that lead to fractional revival can also exhibit perfect state transfer. To 
that end, introduce the Hadamard matrices 


/I 


H = 


V2 


\1 

for IV odd and 

(1 


H = 


V2 


1 \ 


1 1 
1 -1 


- 1 / 


1\ 


1 1 

V2 

1 -1 


\1 

for N even. It is readily seen that 


(6.19a) 


(6.19b) 


e~ iTJ = HU (a) H, 


( 6 . 20 ) 


with a - j - 29 and U(a) the unitary diagonal matrix with elements 


Uij(oi) - Sij 


e ia 

e~ ia 


ij = 

i,j = Lf + 


where [xj is the integer part of x. For M an integer, it thus follows that 


( 6 . 21 ) 


e -iMTj = R u( Ma ) H (6.22) 

since H 2 = 1. Therefore, after a time MT one has |-2#->M(| - 2d). Express now 
the manifestation of fractional revival in the form 
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-iMTJ 


10 } = e‘ 


cos M 




Perfect state transfer will occur if 


(6.24) 


(6.25) 


with Mi an arbitrary odd number since then e lMTJ \0) = e l ^\N) with <f some phase 
factor. Condition (16.25(1 is readily seen to be equivalent to (16. 19(1 : when N is even and 
<5 = 1 — — it is immediate and when N is odd and <5 = 1 + — one uses the properties of 
the cosine to conclude. Hence when (16.2511 is verified these spin chains with FR will also 
exhibit PST at time MT. Take for example the perfectly balanced case of FR which 
occurs when 9 = tt/8, one has then <5 = 1 + — = § for IV odd or <5 = 1 - — = £ for N even 
and it follows that PST will also happen. The evolution will go like this: at time t = T, 
the packet initially at site 0 is revived at 0 and TV, at t = 2 T it is perfectly transferred 
at N, at t = 3 T is is revived again at 0 and N, at t = 4T it perfectly returns to 0, and 
so on. 

Let us now complete our systematic analysis by considering the general case where 
the phase if is arbitrary. As stated at the beginning of the section the polynomials that 
will determine the general Hamiltonians are orthogonal with respect to the weights 
(16.11) associated to the bi-lattices (14.91) . We now understand that we can obtain 
these polynomials in two steps. First, we determine the para-Krawtchouk polynomials 
associated to a bi-lattice with 


<5 


i 4a 

1 ± - 5 

7T 


( 


N odd \ 
N even/ 


(6.26) 


where a is given by 


<7 = 


7 T + rj - £ 
4 


(6.27) 


The Jacobi matrix J is then given by (16.161) and (16.171) with 9 replaced by <x and again 
using the + sign in (16.261) when N is odd and the - sign when N is even. At this point 
we have that e~ lTJ is given by (16.71) or (I6.8[) with 9 again replaced by <r and <f> by 


= f(<v±i)H("?+o. 


( 


N odd 
N even 


(6.28) 


This gives us the polynomials Pi that are associated to the bi-lattice (14.91) but 
are orthogonal with respect to the weights (16.11) with 7 = 1 (corresponding to a 
persymmetric Jacobi matrix) and the of the para-Krawtchouk polynomials. The 
required polynomials that are properly orthogonal against the weights (16.11) are the 
perturbed polynomials Pi defined in (15.361) with 


Co = - 1) 


(6.29) 
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with r a new angle so that 7 - y -1 = - 2 tan 2 r and Ji the recurrence coefficients 
determined in the first step. These polynomials Pi will yield through their recurrence 
relation, the parameters Ji and Bi of the generic chain. The modifications relative to 
the para-Krawtchouk coefficients are given by the formulas (15.13P with 9 replaced by 
r. The determination of e~ lTJ is achieved by conjugating e~ lTJ as given in (16.71) or 
(16. 8 p respectively with the matrix V of (15.7p or (15.6P with 9 replaced by r. Thus are 
determined the Hamiltonians (within the class considered) that have general fractional 
revivals at two sites. Note that as needed, the two-step process has introduced two angles 
a and r. The correspondence with the original parameters 9 and if can be obtained by 
determining explicitly e~ lTJ as indicated before and identifying the coefficients so that 
e -jTj|o) - e %<t> [sin 29 10) + ie 1 ^ cos 2 9 |iV)]. This leads to the following relations 

e^(sin 2 a + i cos 2 a sin 2r) = e ^ sin 29, (6.30a) 

ie 1 ^ cos 2 a cos 2r = cos 29, (6.30b) 


with f given by (16.280 . These conditions are the same for N odd or even provided the 
appropriate (f is chosen. Equation (I6.30bp immediately leads to 


7T 


- <f - if + — + 2nir, 


n e Z. 


The real part of (I6.30ap yields 


(6.31) 


sin 2a - sin 29 sin if, 


(6.32) 


which must be identically satisfied. Upon writing the above equation in the form 

- sin29sinif, (6.33) 


cos 


using (I6.27p . that (I6.32p holds is verified from trigonometric identities having recalled 
the definitions of £ and rj (see the sentence after (14.71) ) and observed that sin£ and sin 7 
must have opposite signs. There then remains from (16.301) the conditions 


cos 2 r = cos 2 #cosec ^ ^ , 

(6.34) 

/ £ — 77 \ 

sin 2 r = sin 29 cos if cosec 1 —-— 1, 

(6.35) 

which determine r. 


7. Conclusion 



Let us summarize our findings. We have completely characterized the XX spin chains 
with nearest neighbor couplings that admit fractional revival at two sites. There 
are two basic ways according to which FR can be realized. One is via isospectral 
deformations of chains with the PST property and the other is by a mirror-symmetric 
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set of couplings corresponding to the recurrence coefficients of the para-Krawtchouk 
polynomials. Hamiltonians with FR at two sites controlled by two arbitrary parameters 
are obtained by compounding these two approaches. The second approach comes with 
a complete set of couplings and magnetic fields while the first approach only sees the 
modification of a few central coefficients of a parent PST chain. The time T for FR 
occurrence doest not depend on the length of the chain. The first method can be 
applied to any PST chain to obtain a chain with FR. Assuming the model is exactly 
solvable to start with, it will remain so under the isospectral deformation. The models 
corresponding to the second way are analytic and may exhibit PST in addition to FR. 
A note is in order here. In principle all chains with FR at two sites can be obtained 
from the generic two-parameter models by surgeries. Indeed, since we are dealing with 
spectra that are finite, any admissible set of eigenvalues can be obtained by removing 
levels from a bi-lattice chosen as large as required. With every such removal, the analytic 
expressions for the chain parameters will become more and more involved thus obscuring 
the exact solvability property. 

It has been indicated that information transfer can be achieved with spin chains 
showing FR at two sites. Knowing that the clone of the initial information will be at 
the end of the chain with definite probability (cos20) 2 at the prescribed time T, the 
end site content at that time can thus be used as input to some quantum process or 
computation with the effect that the final output of the computation will provide the 
right answer with known probability related to (cos2#) 2 . Observe that this probability 
can be tuned by setting correspondingly the chain couplings. Note also that the presence 
of another clone at the site (£ = 0) where the data is entered could be used periodically 
in an experimental or practical context to check that transmission is proceeding without 
alterations since the outcomes at £ - 0 and at £ - N are correlated. It has also been 
pointed out (21 [8, STj that balanced perfect revival can generate entanglement. Indeed it 
is readily observed that for 9 = n/8 the sites 0,1, iV - 1 and N for instance, will support 
at time T the entangled state | t)| i)+| i)| t)- 

Another question has to do with precision. Throughout this paper we have looked 
for situations where FR occurs with probability 1. This could be unduly stringent in 
view of the unavoidable instrumental error for instance. In fact, it would suffice in that 
perspective to consider situations where FR can happen with probability as close to 1 
as desired. This question has been analyzed in the case of full revival and has been 
referred to as almost perfect state transfer (APST) [26] or pretty good state transfer 
[27]. One may assume that the isospectral deformations of a chain with APST will lead 
to chains with almost perfect fractional revival (APFR). Furthermore it has been shown 
in [13] that the para-Krawtchouk chains admit APST for a time T independent of N if 
the bi-lattice parameter 5 is irrational, they should thus admit APFR in those cases too. 
The robustness of FR in the para-Krawtchouk model has also been checked numerically 
in [8]. 

Finally, it would be of great interest to study the possibilities for fractional revival 
at more than two sites. It is known that any unitary matrix can be presented in a form 
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with two diagonals and 2 antidiagonals [2S]; this is related to CMV theory [22] • Assume 
that e~ lTJ is in that form in the register basis |£), where £ - 0,1,..., N. This implies 
revival at up to four sites. A relevant question is to determine the Hamiltonians H 
with their one-excitation restrictions J that will lead to such unitaries. This is likely to 
involve operators beyond the realm of nearest-neighbor interactions. We hope to report 
on this question in the near future. 
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